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Abstract 

We provide a detailed description for power-law scaling FRW cosmological models in Brans-Dicke 
theory dominated by two interacting fluid components during the expansion of the universe. 
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1. INTRODUCTION 

Brans-Dicke (BD) theory is considered as a natural extension of Einstein' s general theory 
of relativity [l| , where the gravitational constant becomes time dependent varying as inverse 
of a time dependent scalar field^which couples to gravity with a coupling parameter u. Many 

9[ can be successfully explained by using this theory. One 



11| for the 



of the cosmological problems 

important property of BD theory is that it gives simple expanding solutions 10| 

scalar field ip(t) and the scale factor a(t) which are compatible with solar system observations 

n n n 

[12]-[l4], in which impose lower bound on uj {\uj\ > 10^) [15|. 

On the other hand, recent observational data give a strong motivation to study gen- 
eral properties of Friedmann- Robert son- Walker (FRW) cosmological models containing more 



than one fluid ( for example, see 16j). Usually the universe is modeled with perfect fluids 
and with mixtures of non interacting perfect fluids. However there are no observational data 
confirming that this is the only possible scenario. This means that we can consider plausible 
cosmological models containing fluids which interact with each other. 



In this work, we follow the authors in 16| but apply the two fluid interaction in FRW BD 
cosmologies. In this case the transfers of energy among these fluids in relation to the BD 
scalar field play an important role in the formalism. There are many cosmological situations 
with the exchange of energy that will be investigated in this work. Although the interaction 
between for example dust-like matter and radiation was previously considered in [l6|, [17 1 
and [18|, the application in BD cosmology gives us new insight to the expansion of the 
universe and the effect of BD scalar field on the subject. 

2. THE MODEL WITH TWO INTERACTING FLUIDS 

We start with the BD action for FRW universes filled with two perfect fluids with energy 
densities pi and p2- The Friedmann equation is given by 



k .ottV ^^^ P1 + P2 



3i7^ + 34 + 3if^-^^ = ^^i:^, (1) 



where k = —1,0, 1. We assume that the two perfect fluids interact through the interaction 
term Q according to 

Pi + 3i7(pi+pi) = g, (2) 

P2 + SH{p2+P2) = -Q. (3) 

Ahhough the nature of the Q term is not clear at all, for Q > there exists a transfer of 
energy from fluid with p2 to the fluid with pi. Also for Q = there are two non-interacting 
fluids, each separately satisfies the standard conservation equation. In the following we 
study the dynamics of the universe for closed, open and flat universes separately. 

2.1. Closed and open power-law interacting cosmologies 

Let us now consider FRW cosmological models with k = —1,1 filled with interacting 
matter sources which satisfy the barotropic equation of state, i.e. 

Pl=llPU P2=72P2, (4) 

where 71 and 72 are equation of state (EoS) parameters. We assume that the scale factor and 
the BD scalar field behave as a{t) oc t" and ip{t) oc a*^ oc t"°, where a and n are constants. 
This implies that H oc at~^. By taking into account the curvature term 3k /a^ of equation 
([1]), we conclude that a = 1 in order to obtain energy density scales over scalar field on the 
RHS of the equation in the same manner as the curvature term in the LHS. Since a oc t has 
no acceleration, the universe will either expand or collapse with constant velocity. 

Now, from equation ([1]) and the resultant equation from the addition of equations ([2]) 
and (j3]) we obtain 

(1 + n + 372) (3(1 + k + n)- n^uj/2) ^^_, 

Pkl{t) = — r 1 , (5) 

3(72 - 7l) 

The interacting term is also found to be, 

g^^) ^ (l+r2 + 37i)(l + n + 372)(3(l + fc + n)-r2V2) ^„_3^ 

3(72 - 71) 
which may be rewritten as 

Q{t) = -(1 + n + 3-fi)Hpki = (1 + n + 3-f2)Hpk2. (8) 



This implies that it is proportional to the expansion rate of the universe and to one of the 
individual densities, so Q ~ t"~'^. 

For the condition ti^uj < 6(1 +n + k), and by defining r"^^ = (3(1 +n + k) — n^(X'/2)t"'~^, 
the densities ^ and (|S]) become, 

(1+72 + 372) 2 /QX 

PkiKV = -— ^r , (9) 

3(72 - 7i) 

(l + n + 371) 2 nm 

^"^"^ = - 3(72-71) " • ^'°^ 

Alternatively for n?uj > 6(1 + n + fc) and taking r"""^ = —(3(1 + n + k) — n^cj/2)t"~^, we 
have 

(l+n + 372) ^^.2 , . 

Pki{r) = ^7- -^r , (11) 

3(72 - 7i) 

(l+n + 371) 2 /.^N 

Pfc2(r) = +^7- -^r . (12) 

3(72 - 7i) 

As can be seen, in figure (la), if we assume that one of the fluids is dust, for different values 

of n, for a particular 71 where pki = Pk2, non of the fluid dominate. Before the crossing point 

Pk2 dominates and after it pki dominates. As n decreases, the domination of pki over pk2 

occurs at larger value of negative 71 while 72 represent dust. Also for those values of n that 

Pk2 < 0, the Weak Energy Condition (WEC) {pki > 0, pk2 > 0) is not satisfied and there 

exist no physical interpretation. In figure (lb) where 72 = 1/3, we have the same argument 

and the only difference is that the curves are shifted smoothly to the right. 

Also, for the condition n'^uj > 6{n + 2) as shown in figures (2a) and (2b), the WEC is not 
satisfied and therefore these is no physical interpretation for this case. 

From WEC and supposing that 72 > 71, and for the condition n'^u < 6(1 + k + n), we 
have 

72>-(n + l)/3, 7i<-(n+l)/3. (13) 

Similarly for the condition n^u > 6(1 + /c + n), we have 

72<-(n + l)/3, 7i>-(n + l)/3. (14) 

From these expressions we conclude that for n > 0, always one of the interacting fluids 
must be either a dark or a phantom fluid. For example the constraints ( !T3|) on the EoS 
parameters imply that Q > 0, so the energy is transferred from a fluid with (—1 > 72 > 





FIG. 1: Plots of pkiT " and pk2T " as functions of 71 with the condition n uj < 6(n + 2) for a) 
72 = 0, b)72 = 1/3 . 

— (n + l)/3) or a phantom (71 < —1) fluid to the dark component whose EoS parameter is 
72 > —{n + l)/3. Also the constraint flT4|) on the EoS parameters imply that Q < 0, so the 
energy is transferred from a matter component whose EoS parameter is 72 > — (n + l)/3 to 
a dark (—1 > 71 > —{n + l)/3) or a phantom (71 < —1) fluid. 

The constant ratio of energies defined by r^ = Pk2/Pki = — Y+n+f^ ^^ another important 
factor that we consider in here. As can be seen, it is a /c independent function of the model 
EoS parameters 71, 72 and n. For cosmological scenarios which satisfy the requirement f ll3p 
, we have r^ = pk2/pki > 1- Thus, pk2 dominates over p^i- On the other hand, for those 
which satisfy (JH]), we have r^ = Pk2/Pki < 1 and pki dominates over pk2- 



2.2. Flat power-law interacting cosmologies 



From equations (JI]), (EJ, and ([3]) for flat FRW model we conclude that the general solution 
for energy distributions are given by. 



_ «(3(1 + n) - n^co/2){in + 3 + 372)^ - 2) ^^^_, 
Pio — ^7 ^ 1 , 

3(72 - 7i) 

a(3(l + n) - n'cu/2){in + 3 + 37i)a - 2) ^„,_, 
P20 — ^7 ^ 1 

3(72 - 7i) 
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FIG. 2: Plots of pkiT " and pk2T " as functions of 71 with the condition n oj > 6(n + 2) for a) 
72 = 0, b)72 = 1/3 . 



The interacting Q-term takes the form 



Q{t) 



a(3(l + ri) - n'^u/2){{n + 3 + 37i)a - 2)((r2 + 3 + 372)0 - 2) 



na—'i 



3(72 - 7i) 

or in terms of the expansion rate of the universe and to one of the individual densities 
((n + 3 + 37i)a-2) ((72 + 3 + 372)^-2) ^ 



(17) 



a 



a 



(18) 



From equation ( IT7j) we conclude that Q ~ t"" ^. For the condition n'^u < 6(1 + n), and by 
definition r""-^ = (3(1 + n) - n2a;/2)r"~2 ^^g gg^g 



, , a(3 + n + 372)a-2^ „^_2 

Pio(r) = —, r r 

3(72 - 7i) 

_ a[(3 + n + 37i)a-2] 2 

P20l7^J — TT ^ ^ 



(19) 
(20) 



3(72 - 7i) 
Alternatively, for the condition n'^u > 6(1 + n), and r'""~^ = —(3(1 + n) — ?T,^a;/2)t'""~^ we 



obtain 



, , a[(3 + n + 372)a-2^ „^_2 

3(72 - 7i) 

, . a[(3 + n + 37i)a-2] 2 
P2o(^) = 57 ^ r 

3 72-71 



(21) 
(22) 



From WEC, for 72 > 71 and the condition r?bj < 6(1 + n), we have 



72 + 3 + 872 
2 



< a < 



n + 3 + 371 
2 



(72 > 7i > 



(71 + 3) 



-00 < a < 



n + 3 + 371 n + 3 + 372 

2 2 

< a < 



, (^ + 3) 
< a < 00, (72 > — > 7i 



n + 3 + 372 n + 3 + 371 

Also for 72 > 7i and the condition r?iij > 6(1 + n), we have 

2 2 



(71 < 72 < 



(n + 3) 



n + 3 + 371 
2 



< a < 



(72 > 7i > 



(71 + 3) 



— cxD, < a < 



71 + 3 + 372 77 + 3 + 371 
2 2 



71 + 3 + 372 ' 

2 , (77 + 3) 

< a < 00, (72 > — > 7i 



< a < 



(71 < 72 < 



(n + 3) 



(23) 
(24) 
(25) 

(26) 
(27) 
(28) 



n + 3 + 371 n + 3 + 372 3 

For n > 0, equations (12 3 p and (1261) are vahd for configurations which include two interacting 
fluids obeying the dominant energy condition (DEC), equations ( l24l) and ( l27l) are valid for 
configurations where one interacting fluid obeys DEC and the other is a phantom fluid, and 
equations (!25l) and (!28|) are valid for the description of two interacting phantom fluids. Now 
let us examine two interesting cases. 

1. Dust-perfect fluid interaction (71 = 0, 72 7^ 0). For the condition n^u < 6(1 + n), and 
from equations. ( TT9l) and ( 120|) . we have 



Pio 



P20 



a((7i + 3 + 372)0 - 2) 



-r 



.na— 2 



372 

-a((n + 3)a-2) ^„^_2 
372 



(29) 
(30) 



Also, for the condition n^w > 6(1 + n), and from equations (12T|) and (!22|) . we have 

pio = _«((^ + 3 + 372)a-2)^„„_2^ ^3^^ 

(32) 



a((n + 3)a-2) 2 
P20 = : 7- 



372 
For the requirement of simultaneous fulfillment of the conditions pio > , P20 > and from 

equations (!23|) - (l25|) the following constraints must be satisfied for the case n^u < 6(n + 2) 



— 00 < a < 



2 
71 + 3 + 372 


< a < 


2 
n + 3' 


2 


2 




^S + 3' 


?7 + 3 + 372 


2 


- < a < 


2 



(72 > > - 



[n + 3). 



(n + 3) 
< a < 00, (72 > -^ — - — - > 0) 



n + 3 + 372 



n + 3 
7 



(0 < 72 < 



[n + 3). 



(33) 
(34) 
(35) 



Also for n^cj > 6(1 + n) and from equations (I26|) - (l28|) we obtain, 



< a < 5 — ^, (72 > > - ^ ^ ' ) (36) 



n + 3 n + 3 + 372 ' 

2 2 , (n + 3) 

-00 < a < , ^ , ^ , ^— < a < 00, 72 > - ' ^ > (37 

n + 3 + 372 n + 3 3 

2 2 , (^ + 3), 

< « < , o , o • (0 < 72 < -^^) (38) 



n + 3 n + 3 + 372 3 

As a specific example, we shall now consider in some detail the dust-radiation interaction 
(71 = 0,72 = 1/3). In this case we have for n'^u < 6(1 + n) 

pio = a[(n + 4)a-2]r""-2, p,=0, 

P20 = -a[in + 3)« - 2]r'^"-2, p^ = ^p^o. (39) 

For the condition n^u > 6(1 + n), we have, 

pio = -a[(r2 + 4)a-2]r""-2, p^ = 0, 

P20 = «[(n + 3)a - 2]r""-^ p2 = ^P2o. (40) 

In figure (3a), where n^u < 6(1 + n), for different values of n, we first have radiation dom- 
inated and then dust dominated. Since the universe in both radiation and dust dominated 
has to be in deceleration phase, we expect that dust domination continues only till a < 1. In 
addition, from figure (3a), n under no condition can be greater than —1.5, since a becomes 
greater than one and the universe starts to accelerate which contradicts the universe decel- 
erating behaviour in radiation dominated and matter dominated era. Further, it can be seen 
that as n decreases, the crossing of the dust and radiation distributions takes place at greater 
values. As shown in figure (3b), there is no physical solution for the case n'^uj > 6(1 + ra). 

2. Phantom-perfect fluid interaction (say 71 = —4/3, 72 7^ 0). For the condition n^o; < 
6(1 + n), we have 

Q;((n + 3 + 372)^-2) 2 

Pio = 5 —. r 

372 + 4 

-a((n-l)a-2) 2 /.^^ 



For the condition n uj > 6(1 + ra 



a((r2 + 3 + 372)0 - 2) _2 



a((n-l)a-2) 2 /.^^ 

8 





FIG. 3: The behaviors of interacting dust distribution, piQ, and radiation distribution, p2o, are 
plotted as functions of a for conditions a) n^uj < 6(1 + n) and b) n?Lo > 6(1 + n). 



As a specific example we shall now consider in some detail the phantom-dust interaction 

(71 = —4/3, 72 = 0). In this case we have for n'^u < 6(1 + n), 



Pio 
P20 



a((n + 3)a-2) 



-r 



.no— 2 



For ri^u > 6(1 + n), 



Pio 
P20 



-a((n-l)a-2) ^^„_2 



a((n + 3)a-2)^„„_2 



(43) 



4 
a((n-l)a-2) 



.no— 2 



(44) 



From figure (4a), one can see that matter domination occurs before phantom domination. 
Further, it can be seen that as n decreases, the crossing of the phantom and dust distributions 
takes place at greater values of a. Moreover, for a = 1 when the universe starts accelerating 
and phantom dominates over matter, the value of n is one. Figure (4b) shows no physical 
interpretation for the case n^u > 6(1 + n). 
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FIG. 4: The behaviors of interacting phantom distribution, pio, and matter distribution, p20) for 
conditions a) rt^uj < 6(1 + n) and b) n^LO > 6(1 + n). 

3. THE EFFECTIVE FLUID INTERACTION 



In this section we study the conditions under which the two interacting sources are 
equivalent to an effective fluid described as the interaction of two-fluid mixture. This can 
be made by introducing an effective pressure Peff, i-e. 



Pl+P2 = llpl + 72P2 = Peff, 



which has an equation of state given by 



Peff = lPeff = 7(Pl+P2), 



(45) 



(46) 



where 7 is the effective EoS parameter. Note that the equation of state of associated effective 
fluid is not produced by physical particles and their interaction [17]. Making some algebraic 
manipulations with equations (H5l) and ( H6l) we find that the the effective EoS parameter, 
7, is related to the parameter a by, 

2 + n 



ri^oj + n^ — n 



7 



+ 



3(n + 1) - n2a;/2 a[?,{n + 1) - n'^uj/2]' 
We also find the behavior of the BD parameter u in terms of n (n 7^ 0)in: 
Inflation era [a » 1, 7 = —1): u ■ 



(47) 



-1 + i 

n 
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Radiation dominated era {a = 1/2, 7 = 1/3): u = —3 — - 

Matter dominated era {a = 2/3, 7 = 0): u = —2 — - 

Dark energy dominated era {a ^ 1.1, 7 = —1): cu = —I + ^ + -^ 

It shows that in the inflation and cosmological constant era the behaviour of u in terms 

of n is different from in radiation and matter dominated era. In case of n = { or (p = (J)q) 

eq. p7|) becomes, 

2 — 3a ,^. 

7 = ^r— , 48 

6a 

in compatible with the standard cosmology. That is, in radiation era where a = 1/2, the 
EoS parameter becomes 1/3, in matter dominated era where a = 2/3, the EoS parameter 
becomes zero, in inflation era where a — )■ 00, the EoS parameter becomes —1 and in the 
current dark energy era where a is just greater than one, the EoS parameter becomes about 
— 1/3 which is compatible with the observational data. 

4. DISCUSSION 

In this paper we have provided a detailed description for power law scaling cosmological 
models in the case of a FRW universe in BD theory dominated by two interacting perfect fluid 
components during the expansion. We have shown that in this mathematical description it 
is possible for each fluid component to require that the WEC may be simultaneously fulflUed 
in order to have reasonable physical values of EoS parameters in the formulation. So from 
the required conditions we may gain some insights for understanding essential features of 
two fluid interactions in power law FRW BD cosmologies. In closed or open universes, if 
one of the fluids is dust or radiation, the second fluid has to be dark energy or phantom. 
As the power of BD scalar fleld, n, a function of time, decreases, the second fluid changes 
its behaviour from dark energy to phantom. If one of the fluids is radiation, the change of 
behaviour of the second fluid, as n decreases, is slower in comparison to the dust fluid. 

In flat universes, if the two fluids are dust and radiation, their behaviuor as a function of 
scale factor power of time, a, changes with n. It has been shown that as n decreases, the 
domination of dust over radiation occurs for smaller values of a. It also gives us a constraint 
on n not to be greater than values that makes the universe accelerate. Similarly, if the 
two fluids are phantom and dust, again their behaviour as a function of a changes with n. 
Firstly, for a particular n, the domination of dust occurs before the domination of phantom. 

11 



Secondly, as n decreases the domination of phantom occurs at greater values of a. Thirdly, 
n is constrained such that in dust dominated era a has to be less than one for a decelerated 
universe. 

Finally, we investigate the effective fluid interaction in the model. We find an expression 
for the effective EoS parameter, 7, as a function of a, n and u. The behaviour of the BD 
parameter, w, as a function of n in different epoch of the universe expansion is given. It 
has been found that uj is inversely (negatively) proportional to n in radiation and matter 
dominated era, whereas inversely (positively) proportional to n in inflation and dark energy 
dominated era. For the constant scalar field the effective EoS parameter turns out to be the 
one expected in standard cosmology. 
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